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ABSTRACT 

A design optimization method for turbopumps of 
cryogenic rocket engines has been developed. 
Multiobjective Evolutionary Algorithm (MOEA) is 
used for multiobjective pump design optimizations. 
Performances of design candidates are evaluated by 
using the meanline pump How modeling method based 
on the Euler turbine equation coupled with empirical 
correlations for rotor efficiency. 

To demonstrate the feasibility of the present approach, 
a single stage centrifugal pump design and multistage 
pump design optimizations are presented. In both 
cases, the present method obtains very reasonable 
Pareto-optimal solutions that include some designs 
outperforming the original design in total head while 
reducing input power by 1%. Detailed observation of 
the design results also reveals some important design 
criteria for turbopumps in cryogenic rocket engines. 
These results demonstrate the feasibility of the EA- 
based design optimization method in this field. 

INTRODUCTION 

While the budget for space development programs has 
drastically shrunk in most countries, recent and future 
space missions increasingly demand high performance 
and reliable rocket engine systems and components, 
such as turbopumps. Progress in computational fluid 
dynamics (CFD) methods and development of 
powerful computational facilities have contributed to 
the reduction in required cost and time to develop 
advanced turbopump designs. The design process still 
largely depends on experienced designers. Therefore, 
numerical design methods coupled with CFD, which 
are capable of efficiently developing advanced 
turbopump designs, can greatly reduce such 
dependency. 


Among numerical optimization algorithms, gradient- 
based methods are long-standing and most widely used 
approaches 1 ' 3 . These methods use the gradient of an 
objective function with respect to changes in design 
variables to calculate a search direction using steepest 
descent, conjugate gradient, quasi Newton techniques, 
or adjoint formulations. The solution obtained by these 
methods will be a global optimum, only if the objective 
and constraints are differentiable and convex . 
Unfortunately, the distribution of an objective function 
of real-world design problems is usually multimodal 
and one could only hope for a local optimum 
neighboring the initial design point. Therefore, to 
determine the global optimum, one must optimize from 
a number of initial points and check for consistency in 
the optima obtained. In this sense, the gradient-based 
methods are not robust. 

Evolutionary Algorithms (EAs, for example, see [5]) 
are emerging design optimization algorithms modeled 
on the mechanism of natural evolution. EAs search 
from multiple points, instead of moving from a single 
point. In addition, they require no derivatives or 
gradients of the objective function. These features lead 
to robustness and simplicity in coupling with any 
evaluation codes. Parallel efficiency also becomes very 
high by using a simple master-slave concept for 
function evaluations, if such evaluations consume most 
of CPU time. Design optimization using CFD is a 
typical case. Application of EAs to multiobjective 
design problems is also straightforward because EAs 
maintain a population of design candidates in parallel. 
Due to these advantages, EAs are unique and attractive 
approach to real-world design optimization problems. 
Recently, EAs have been successfully applied to 
aerospace design optimization problems' ' . 
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The objective of the present study is to develop and 
demonstrate a design optimization method for 
turbopumps used in the cryogenic rocket engines. The 
Multiobjective Evolutionary Algorithm (MOEA) will 
be used for multiobjective optimization of pump 
designs. Performances of design candidates will be 
evaluated by using the meanline pump flow modeling 
method based on the Euler turbine equation coupled 
with empirical correlations for rotor efficiency. Present 
approach will be applied to centrifugal and multistage 
turbopump design optimization problems. 

EVOLUTIONARY ALGORITHMS 

EAs mimic mechanism of natural evolution, where a 
biological population evolves over generations to adapt 
to an environment by selection according to fitness, 
recombination and mutation of genes (Fig.l). When 
EAs are applied to optimization problems, individual, 
fitness, and genes usually correspond to a design 
candidate, an objective function value, and design 
variables, respectively. One of the key features of EAs 
is that it searches from multiple points in the design 
space, instead of moving from a single point like 
gradient-based methods do. Furthermore, these 
methods work on function evaluations alone and do not 
require derivatives or gradients of the objective 
function. These features lead to the following 
advantages: 

1) Robustness: EAs have capability of finding a 
global optimum, because they don't use function 
gradients that direct the search toward a local 
optimum. In addition, EAs have capability to 
handle any design problems that may involve non- 
differentiable objective function and/or a mix of 
continuous, discrete, and integer design 
parameters. 

2) Capability of sampling various Pareto-optimal 

solutions in parallel: Real-world design 

optimization problems typically involve multiple 
and often competing objectives. The solution to 
such problem is not a unique optimal solution, but 
a set of compromised solutions, largely known as 
Pareto-optimal solutions. Each of those solutions is 
optimal in the sense that no improvement can be 
achieved in one objective component that does not 
lead to degradation in at least one of the remaining 
components. Therefore, primary goal of a 
multiobjective optimization problem is, unlike that 
of a single objective optimization, to find various 
Pareto-optimal solutions to show the precise 
tradeoff information among the completing 
objectives. By maintaining a population of 
solutions and introducing the concept of Pareto- 
optimality, EAs can uniformly sample various 
Pareto-optimal solutions in parallel. 


3) Suitability to parallel computing: Since EAs are 
population-based search algorithms, all design 
candidates in each generation can be evaluated in 
parallel by using the simple master-slave concept. 
Parallel efficiency is also very high, if objective 
function evaluations consume most of CPU time. 

4) Simplicity in coupling evaluation codes: As these 
methods use only objective function values of 
design candidates, EAs do not need substantial 
modification or sophisticated interface to 
evaluation codes. If an all-out re-coding were 
required to every optimization problem like the 
adjoint methods, extensive validation of the new 
code would be necessary every time. EAs can save 
such troubles. 

The present MOEA uses floating-point representation, 
where an individual is characterized by a vector of real 
numbers. It is natural to use the floating-point 
representation for real parameter optimization 
problems instead of binary representation because it is 
conceptually closest to the real design space, and 
moreover, the string length is reduced to the number of 
design variables. Fonseca's Pareto-based ranking 
method 50 is used for fitness assignment where an 
individual's rank corresponds to the number of 
individuals in the current population that are better than 
the corresponding individual in every objective 
function. To maintain diversity in the population, a 
standard sharing function 10 is incorporated. As the 
elitism, the best-N selection 1 1 is incorporated, where 
the best N individuals are selected for the next 
generation among N parents and N children based on 
Pareto-optimal ity“ so that the Pareto-optimal solutions 
will be kept once they are formed. Parents are selected 
from the best N individuals randomly. To generate new^ 
design candidates from these parents, the blended 
crossover (BLX-a) is used, which is the most common 

approach for recombination of two parents represented 
by a vector of real numbers proposed by Eshelman and 
Schaffer 1 ". In this approach, children are generated on a 
segment defined by two parents, but the segment may 
be extended equally on both sides determined by a user 
specified parameter a. Thus, a child solution is 
expressed as: 

Child I = 7 • Parent) + ( 1 — y) • Parent 2 ( 1 ) 

Child 2 = (1 — y) • Parent) + y ■ Parent 2 (2) 

w r here 

y= ( l + 2a) u - a ( 3 ) 

Child 7, Child2 and Parentl , Parentl denote design 
parameters of the children and parents, respectively, u 
is uniform random number in [0,1]. Schematic view of 
BLX-a is shown in Fig. 2. When an EA is applied to a 

design optimization problem, what is important is 
balance of two conflicting goals: exploiting good 
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solutions and exploring the search space 13 . Thus, BLX- 
0.5 is used in which both exploration and exploitation 
are carried out equally. Since the strong elitism is used, 
high mutation rate of 0.2 is applied and a random 
disturbance is added to the parameter in the amount up 
to i 20% of the design space. Population size and 
maximum number of generations are set to 100 and 90 
for the centrifugal pump design, 100 and 120 tor the 
multistage pump design. Unbiased initial population is 
generated by randomly spreading solutions over the 
entire design space in consideration. 


PUMP PERFORMANCE EVALUATION 
Total head and required input power ot pump design 
candidates are evaluated by using a one-dimensional 
meanlinc pump flow modeling method 14 , which 
provides a fast capability for modeling turbopumps 
within rocket engines. The components of the inlet and 
exit fluid velocity triangles are calculated at the hub, 
mean and tip locations along the rotor blades. The 
meridional velocity of the fluid at the rotor leading 
edge root-mean-square diameter C\u [ft/sec] is defined 
by equation (4). 


Cm i 


1 44/1? 

pA 


( 4 ) 


where 

m\ Mass flow [lbs/sec] 

p/; Fluid density at leading edge [lbs/ft'] 

A/: Flow area at leading edge [inch 2 ] 

Flow area is calculated from the input flow path 
dimensions. 

A = ApB(R huh + R, p )- block ] ( 5 ) 


where 

A: Boundary layer blockage factor 

B : Blade span from hub to tip [inch] 

R, wh : Radial distance from pump centerline 
at hub [inch] 

R ljp : Radial distance from pump centerline 

at tip [inch] 

The metal blockage of the rotor block is calculated by 
eq. (6). 


block = 


thk-BZ 
sin /3 


( 6 ) 


where 

thk: Normal blade thickness [inch] 

Z: Blade number 

ft: Relative angle from tangential [degree] 

The tangential component of velocity entering the rotor 
is calculated in terms of the swirl angle of the flow a } 


by equation (7). 

Cm - C m / tan(a, ) (7) 


The meridional and tangential components of absolute 
fluid velocity at the rotor trailing edge are calculated by 
equations (8) and (9). 


C M2 — 


144 m 
P 2 A 2 


( 8 ) 


Cu2 = C :+ Wy 2 


( 9 ) 


where 

p : : Fluid density at trailing edge [lbs/ff ] 

Ay. Flow area at trailing edge [inch 2 ] 

U< Blade tangential velocity at trailing edge 
[ft/ sec] 

Tangential component of relative fluid velocity 


at trailing edge [ft/sec] 

Flow area at trailing edge is calculated by eq. (5). The 
blade tangential velocity U and tangential component 
of the fluid relative velocity M V: are given by equations 
(10) and (11), respectively. 

U= 2kR -— (10) 

770 


Wv2= Cw2-tanj8 2 +t/_’f7-<7? (11) 


where 

R: Radial distance from pump centerline [inch] 

A r : Shaft rotative speed [rpm] 

/T: Relative angle from tangential at trailing edge 
[degree] 

The slip factor a is defined by 


G ~ 


slip 

U 2 


( 12 ) 


The slip is the difference between the theoretical and 
absolute fluid tangential velocities. For centrifugal 
impellers, Pfleiderer correlation to geometry 1 '' is used 
to calculate the slip factor a. A default slip factor of 


0.95 is used for inducers. 

The head rise through the rotor is calculated iteratively 
from the Euler turbine equation coupled with empirical 
correlations for rotor efficiency 

ft ,, 


H, =. 


£ 


where 

Hy Head rise through the rotor [ft] 

rj hy(i : Rotor hydraulic efficiency 

g c : Gravitational constant, 32.l74[lhm-ft/lbf-sec‘] 

The rotor hydraulic efficiency is obtained^ from 
empirical correlations to rotor-specific speed 1 . The 
total pressure and static pressure at the rotor exit are 
estimated from the rotor head rise by equations (14) 
and (15). 


Pc 


Hyp t _ 2 
1 44 


+ Pa 


(14) 
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Ps2~ P,2- 


C;-p 2 
2- 144- g r 


where 


( 15 ) 


P,j\ Total pressure at the leading edge [psia] 

P ,: : Total pressure at the trailing edge [psia] 
pi < Average density of the fluid from the leading 

edge to the trailing edge [lbs fr ] 

P s < Statie pressure at the trailing edge [psia] 

Fluid absolute velocity at trailing edge C 2 [ft/sec] is 
defined by 


c = 



(16) 


Total pressure at the discharge of the last stage P t4 
[psia] is given by the following equation 

P ,4 =P l 2-° ) 2-4 ( P ,2~ P ,2) ( 17 ) 

where the design point total pressure loss coefficient of 
the diffusion system is assumed to be known and is 
input in terms of a normalized loss coefficient (0 24 . The 


total head rise through pump is calculated by 

„ _ 144 ( 7^ 4 - P n ) 

H 4 ( 18 ) 


where 

H 4 \ Total head rise through pumpfft] 

P , 4 : Total pressure at the pump exit [psia] 
p /. 4 : Average density of the fluid from the inlet 
to the discharge [lbs ft '] 


The input power required to drive the rotor is 
calculated from the head rise through the rotor, mass 
flow, rotor hydraulic efficiency, mechanical efficiency, 
volumetric efficiency and disk pumping loss as 


( 


input = 


m • H 2 

-77™ 


-+PL 


iisk 


^7 /»{*( h 


(19) 


where 


input: Input power [hp] 

7] Mechanical efficiency 


Vj voi : Volumetric efficiency 
PLj lsk \ Disk Pumping loss [hp] 


The mechanical efficiency is assumed to be 0.98 and 
the volumetric efficiency is based on internal leakages 
and is expressed as the ratio of leakage to the inlet 
flow. The disk pumping loss is calculated from 
empirical correlations to geometry, fluid density at 
rotor trailing edge, and the shaft rotative speed 16 . 
During the calculation, local static pressure at the rotor 
tip is compared to the local vapor pressure to check for 
the cavitation inception point. 

To estimate off-design total head and required input 
power, the empirically derived variation of slip factor 
and rotor efficiency as a function of flow'-speed ratio F 
is used. Correction factor is also applied to the total 
pressure loss coefficient of the diffusion system as a 
function of loading parameter L. 


— — — = 1.534988-0.6681688- F + 0.077472 F 2 + 0.0571508-F' 

a, 

design 

— — — = 0.86387 + 0.3096 -F-0.1 4086 • F 2 - 0.029265 • F ? 

^1 /mi .design 

— — — = 1.8151-1. 83527 -Z. + 0. 8798 ■ Zr +0.1 8765 -I 3 

® design 


( 20 ) 


(21) 


( 22 ) 


The loading parameter is defined in terms of the velocities 
at the vaneless diffuser exit and the velocity at the 
diffusion system throat. 


L = 


C. 


iham 


VC? + Tw? 


(23) 


where 

C lhtoaI : Fluid absolute velocity at the diffusion system 
throat [ft sec] 

(\ p Tangential component of fluid absolute velocity 
at vaneless diffuser exit [ft/sec] 

Cvm: Meridional component of fluid absolute 
velocity at vaneless diffuser exit [ft/sec] 


The velocity at the diffusion system throat is defined by 
the equations (24) and (25). 

144 a ?? 


Cy/ 3“ 

P3^3 


(24) 


c ( „ =(1-CJ 2 „3)C L: 


R2 
R 3 


(25) 


where the pressure loss coefficient at the diffuser exit ok* 


is assumed to be 0.1. Fluid velocity at the throat is given 
by the equation (26). 


C throat " 


144 aa ? 

P } A \hnm 


(26) 
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CENTRIFUGAL PUMP DESIGN 

First, redesign of a single-stage centrifugal pump, M-l 
oxygen scaled tester is demonstrated. Objectives of the 
present design problem are maximization of total head and 
minimization of input power at a design point. These 
objectives are competing and therefore the solution to this 
optimization problem is Pareto-optimal solutions. 

The design point is shaft rotative speed ot 5,416.7[rpm], 
total temperature of the fluid entering the pump of 545 
[Rankine], total pressure of the fluid entering the pump of 
50.0 [psia]. and mass flow into the pump of 188.7 
[lbm/sec]. 

Design parameters are rotor leading edge tip radius ( R lip i ), 
rotor trailing edge radius (R:), volute tongue radius (/?A 
blade span at trailing edge (£_>), blade span at volute 
tongue (£<), axial length of the blade at the root-mean- 
square diameter ( S ), number of blades (Z„), blade thickness 
(thk\ blade trailing edge angle at the hub. root-mean- 
square radius, and tip ($,„/, , /3 , Pn r ) as shown in Fig. 3. 

Table 1 presents present design spaces. 

Total head and input power of Pareto-optimal designs, 
original design, and all other design candidates are 
illustrated in Fig. 4. Designs that have cavitation are 
eliminated from the figure. Present Pareto-optimal 
solutions successfully displays tradeoff information 
between maximization of the total head and minimization 
of the input power. Such tradeoff information is very 
helpful to a higher-level decision-maker in selecting a 
design with other considerations. Among the Pareto- 
optimal solutions, some designs outperform the original 
design in the total head while reducing the input power by 
1 %. 

Figure 5 compares overall performance maps of the 
original design, the highest total head design, the lowest 
input power design, and a compromised design that 
overcomes the original design in both objectives. 
Moreover, the compromised design is seen to improve the 
exit pressure in all off-design conditions. The design 
parameters of these designs are shown in Table 2. 

The absolute flow velocity at rotor exit hub is shown in 
Fig. 6. This figure indicates the optimum designs have 
small exit flow velocity, which contributes to minimizing 
the total pressure loss in the diffusion system. By 
minimizing the total pressure loss in the diffusion system, 
designs can improve their total head rise. To minimize the 
exit flow velocity, the optimum designs have small slip 
factor values, i.e., large slip than others. Actually, all three 
optimum designs in Table 2 have small axial length, while 
the high head design and the compromised design 
maximize R : to increase slip due to the inertial effect (Low 
input design minimize Rj to minimize input power). The 
low input design and the compromised design also 
minimize their blade angle at hub and tip to reduce 
absolute fluid velocity at rotor exit (the high head design 
maximize blade trailing edge angles to improve its total 
head). However, it is known that non-uniform radial 


velocity around the periphery of the impeller due to large 
slip degrades its head rise. Therefore, this effect should be 
counted in the future study. 

Figure 7 shows the total pressure loss coefficient ot the 
diffusion system of the designs. Pareto-optimal solutions 
successfully minimize it to increase their total head rise. 
According to these detailed observations of the results, 
present MOEA obtained reasonable Pareto-optimal 
solutions, which ensure the feasibility of the present design 
optimization approach in rocket engine pump designs. 

MULTISTAGE PUMP DESIGN 

Next, present design optimization method is applied to the 
redesign of RL10A-3-3A liquid oxygen pump consisting 
of one inducer and a single centrifugal impeller, followed 
by a vaneless diffuser and conical exit volute. The 
objectives are maximization of total head and 
minimization of input power at the design point, which is 
shaft rotative speed of 12,900[rpm], total temperature of 
the fluid entering the pump of 175 [Rankine], total 
pressure of the fluid entering the pump of 40.0 [psia], and 
mass flow into the pump of 40.0 [lbm/sec]. Design 
parameters and the corresponding parameter ranges are 
shown in Fig. 8 and Table 3, respectively. 

Figure 9 shows total head and input power of Pareto- 
optimal designs, original design, and all other design 
candidates that have no cavitation. Though this design 
optimization problem involves two stages and a large 
number of design parameters, the present MOEA finds 
reasonable Pareto-optimal solutions including some 
designs that improve both total head and input power by as 
much as 1%. 

Figure 10 shows overall performance maps of the original 
design and optimized designs. The compromised design 
improves the exit pressure in all off-design conditions. The 
design parameters of these designs arc show n in Table 4. 
Figures 1 1 and 12 illustrate the head rise and the required 
input power of the first and the second stages. Because the 
exit of the first stage connects with the inlet of the second 
stage directly, the relation between the head rise and input 
power becomes linear. The optimum designs increase their 
head rise through the first stage because the slope of the 
curve consisting of Pareto-optimal solutions in Fig. 12 is 
steeper than that of the line in Fig. 1 1 . To increase the head 
rise through the first stage, the highest head design and the 
compromised design increase R np , p:hub> p 2 rms* an d p:n P - 
Another interesting thing is that the second stages of the 
Pareto-optimal designs are not optimal by themselves, 
especially in the high head region. This is due to the 
nonlinear interactions between the first and the second 
stages, which make a multistage pump design very 
difficult. 

Figure 13 shows fluid velocity at rotor exit hub. The 
Pareto-optimal designs have a small fluid velocity at exit 
like the single stage pump design. However, the optimum 
designs that have the total head of 950-1200 inches do not 
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minimize their total head. This is probably due to 
complicated interactions between the first and second 
stages. 

Figure 14 is the total pressure loss coefficient of the 
designs. This figure is also interesting because the optimal 
designs in low total head region minimize their total 
pressure loss coefficient but the optimal designs in the high 
total head region maximize it. More work is necessary to 
understand the multistage pump designs. 

SUMMARY 

In the present study, a design optimization method for 
cryogenic rocket engine turbopumps has been developed. 
Multiobjective Evolutionary Algorithm is used for the 
multiobjective optimization of pump designs. 
Performances of design candidates are evaluated by using 
the meanline pump flow modeling method, which is based 
on the Euler turbine equation coupled with empirical 
correlations for rotor efficiency. 

To demonstrate the feasibility of the present approach, 
single stage centrifugal pump design and multistage pump 
design optimizations are presented. In both cases, present 
method obtains reasonable Pareto-optimal solutions that 
include designs outperforming the original design in total 
head as well as input power by 1%. Detailed observation 
of the design results also reveals some important design 
policies in turbopump design of cryogenic rocket engines. 
These results ensure the feasibility of EA-based design 
optimization method in this field. 
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Figure 1. Flowchart of the present evolutionary 
algorithm. 
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Figure 3. Design parameters of the centrifugal pump 
design problem. 


Table 1. Design parameter ranges of the centrifugal 
pump design problem. 
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B, 

S 


[inch] 

[inch] 
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Figure 5. Pump overall performance map. 
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Table 2. Pareto-optimal designs of the centrifugal 
pump design problem. 


design variables 

&ifip 

[inch] 

r 2 

[inch] 

lineh] 

b 2 

[inch] 

B 

[inch] 

S 

[inch] 

High head design 

3.50 

5.60 

6.05 

0.846 

0.871 

3.74 

Low input design 

3.71 

5.00 

6.08 

0.701 

0.859 

3.74 

Compromised design 

3.92 

5.59 

5.63 

0.730 

0.878 

3.76 

Original design 

3.66 

5.34 

5.91 

0.814 

0.908 

4.00 


design variables 

P hub 

[degs.] 

[degs.) 

(degs.) 

thk 

[inch] 

Z* 

High head design 

44.9 

36.0 

44.6 

0.100 

18 

Low input design 

25.3 

33.2 

26.1 

0.065 

6 

Compromised design 

26.7 

37.2 

26.4 

0.042 

7 

Original design 

35.0 

35.0 

35.0 

0.050 

12 



0 I l _ t _£i i I 

600 700 800 900 1000110012001300 
Head [inch] 


Figure 6. Flow velocity at second stage exit hub. 



head [inch] 

Figure 7. Total pressure loss coefficient of the 
centrifugal pump designs. 



Figure 8* Design parameters of the multistage pump 
design problem. 

Table 3. Design parameter ranges of the multistage 


pump design problem. 
<lst stage> 


design variables 

[inch] 

8 2huh 

[inch] 

ft 2huh 

[degs.] 

P jhns 
[degs,] 

P Jv- 
[degs] 

thk : 

[inch] 

Upper boundary 

1.1H 

0.53 

50.0 

40.0 

35.0 

0.08 

Loner boundan ■ 

1.08 

0.43 

35.0 

25.0 

20 0 

0 03 


<2nd staoe> 


design variables 

R. 

[inch] 

[inch] 

B. 

[inch] 

By 

[inch] 

thk , 

[inch] 

S 

[inch] 


p, 

[degs.] 

l ppei boundary 

2.20 

0.15 

0.35 

0.50 

0.08 

1.00 

16 

90.0 

Lower boundary 

2.00 

0.05 

0.15 

0.30 

0.03 

0.80 

8 

60.0 



800 900 1000 11001200 1300 1400 1500 
Head [inch] 


Figure 9. Objective function values of the multistage 
pump designs. 
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Figure 10. Pump overall performance map. 
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Figure 11. First stage performances of the multistage 
pump designs. 


Table 4. Pareto-optimal designs of the multistage 


pump design problem. 
<lst stage > 


design variables 

[inch] 

R.'kub 

[inch] 

fi Jmb 

[degs] 

[degs] 

/»* 
|degs ] 

thk : 

[inch] 

High head design 

1.17 

0.517 

49.9 

3X.5 

34.9 

0,07 XX 

Low input design 

1.17 

0,525 

42.9 

37.3 

2X.0 

0.037ft 

C o mpromised design 

1.17 

0.514 

4X.3 

40.0 

33.4 

0.0520 

Original design 

1.13 

0.4X0 

43.0 

27.3 

21.6 

0.0400 


<2nd staae> 


design variables 

[inch] 

(R 4 -fid 

[inch] 

B s 

[inch] 

b 4 

[inch] 

thk } 

[inch] 

S 

[inch] 

Z„, 

p i 

[degs.] 

High head design 

2.20 

(L 1350 

0.340 

0.301 

0.0311 

0.X55 

16 

X7.4 

Low input design 

2 00 

0.054 X 

0.152 

0.393 

0.0791 1 

0.X04 

X 

64.4 

Compromised design 

2.16 

0.060X 

0.333 

0.370 

0 0439 

0X12 

9 

77.X 

Original design 

2.10 

0.0950 

0.251 

0.400 

0(hm 

0.X7X 

12 

90.0 


Oh 

33 
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Figure 12. Second stage performances of the 
multistage pump designs. 
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Figure 13. Flow velocity at second stage exit hub. Figure 14. Total pressure loss coefficient of the 

multistage pump designs. 
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